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Abstract 

By the use of Perron-Frobenius theory, simple proofs are given of the Fundamental 
Theorem of Demography and of a theorem of Cushing and Yicang on the net repro- 
ductive rate occurring in matrix models of population dynamics. The latter result is 
further refined with some additional nonnegative matrix theory. When the fertility 
matrix is scaled by the net reproductive rate, the growth rate of the model is 1. More 
generally, we show how to achieve a given growth rate for the model by scaling the 
fertility matrix. Demographic interpretations of the results are given. 



1 Introduction 

A standard matrix model of population dynamics is given by a sequence of nonnegative 
vectors xq, x±, . . . of fixed length n defined by 



for a given nonzero xq where P is an n x n matrix with nonnegative entries. As usual, 
we assume that 



where T and F are nonnegative nonzero matrices such that all the column sums of T 
are not larger than one: 



x k = Pxjc-i, A; = 1,2,..., (1) equl 



P = T + F (2) equ2 



£»=!,... ,„ Uj < 1, j = 1,... ,n. (3) equ3 
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The rath entry of the vector represents the population in the mth class at the time 
epoch k. The matrix P = T + F is known as the projection matrix; the matrices T 
and F are known as the transition matrix and the fertility matrix resp., so that the 
entry of T represents the fraction of the individuals in the jth class that will 
survive and move to the ith. class in a unit interval of time, and the entry of F 
represents the number of newborns in the ith class that descend from one individual 
in the jib. class in a unit interval. Because of these demographic interpretations, one 
sees why the column sums of T are always less than or equal to 
of this mat rix mo del for Leslie matrices, introduced by Bernardelli 
and Leslie IfLesl], has a long history R. We concentrate on properties of the general 



The specia' 
~ern|, Lewis 



case 
an- 



ew 



case introduced above. 

The following argument shows that a further assumption is warranted for our matrix 
model. The spectral radius p(T) of a matrix T is the maximum of the moduli of 
its eigenvalues. Since we wish to exclude the possibility of an immortal population, 
we shall always assume in the rest of the paper that lim^oo T k XQ = for all initial 
populations (nonnegative vectors) xq. By examining the effect of T k on the standard 
basis of unit vectors, this condition may be shown to be equivalent to lim^oo T k = 0. 
In turn, it is known that this condition is equivalent to 



P{T) < 1, 

Theorem 3.5], ||J], Theorem J3. 6 L2]. We shall thi 



(4) 



( . •■'••]■ fev'u'l •"i«"u2" J ..^ always assume that our model 

satisfies (3j) together with (jjjjand (Qj; the assumption (j§|) is not usedjn r o_ur proofs. We 
remark that (yj implies the weaker condition p(T) < 1, but that (Qj does not imply 
(Hj as is easily shown by examples. 



Under the assumption 



prt 

we have 



e.g 



mort 



(I -T)- 1 = I + T + T 2 + ... , (5) |defq 

.g. fiw|, Theorem 6.1], |i| Corollary 5.6.16]. Let Q = F(I - T)' 1 . Then we have 

Qx = Fx + FTx + FT 2 x + . . . , 



which represents the distribution with respect to state-at-birth of all newborn descen- 
dants accumulated during the entire lifespan of the population Xq, see | Cush , p.&land 



p. 71], where the matrix Q is called the "next generation matrix" . Following |C as 2 



|!Hy 126], w e, call t he spectral radius p(Q) of Q the net reproductive rate of the model, see 
I CY| and I flCusb , p. 8], for a definition in this generalit y. As usual it is denoted by Rq. 
Other names for Rq are "net reproductive number" ( f|Cusb|l ) and "basic reproduction 
ratio" (pip. 

The purpose of this note is to explore some applications of the Perron-Frobenius theory 
to results in Population Dynamics; particularly we wish to study the role of the net 
reproductive rate. In Section 2 we state needed standard results from this theory 
which will be used throughout the paper. We first apply the theory to obtain the 
Fundamental Theorem of Demography for a mo del, w ith a pr^nfitjye projection matrix. 



If], or jClg 



This approach is well-known, see for example 
jn:sl[ acknowledges the contribution of Ibcxi! and |^s2| mentions 



j 0as2 1 . In Section 3, we 



|bew ] , but l |jcrn| and fiL°cw| do not 



refer to each other, which is extraordinary as both authors were at University College, Rangoon, Burma, in 
1940 - 41 when their papers were submitted. 
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consider the net reproductive rate. We state and give a short proof of a somewhat 
stronger form of an interesting comparison theorem on the net reproductive rate due 



to Cushing and Yicang, see |CY, Theorem 3] and Ij Cush , Theorem 1.1.3]. The theorem 
is further refined in Section 4 using a result from graph theoretic Perron-Frobenius 
theory. For a model with an irreducible projection matrix, the net reproductive rate 
may be viewed as a factor producing a model with growth rate 1 when one scales only 
the fertility matrix by this constant factor. More generally, given a positive s subject 
to one restriction, we determine a q(s) as a function of s such that the growth rate of 
the model is s when the fertility matrix is scaled by q(s). In our last section, we give 
demographic interpretations of our results. 

2 Perron-Frobenius Theory and the Fundamen- 
tal Theorem of Population Demography 

In this sect ion w e give the sketch of ^gro of of the fundamental theorem of demog- 



raphy, see IflCush, Theorem 1.1.2] or |Cas2, p. 86] using the Perron-Frobenius theory 



of nonnegative matrices. This theory is also needed for applications in subsequent 
sections. 

A nonnegative matrix is irreducible if it is not the lxl zero matrix and it is not 
permutationally similar to a matrix of the form 

A 22 

for nontrivial square matrices An and A 22 . The celebrated Perron-Frobenius theorem 
for irreducible matrices can be stated as follows. 

Theorem 2.1 Let P be an irreducible nonnegative matrix. Then 

(a) The spectral radius p(P) of P is positive and it is an algebraically simple eigenvalue 
of P with corresponding left and right positive eigenvector, which are unique up to scalar 
multiples. 

(b) The spectral radius of P is the unique eigenvalue with a nonnegative eigenvector. 

(c) The spectral radius of the matrix P increases (strictly), resp. decreases, if any entry 
of it increases, resp. decreases. □ 

One may see |w|, Theorem 10.7], |hj|, Theorem 8.4.4, Problem 15, p. 515], |v|, Theorem 
2.1], |G], Theorem 2, p. 53] for proofs. 

Following a common practice in matrix literature, we call the spectral radius p{P) of a 
nonnegative matrix P the Perron root of P. In some fields of mathematics this term is 
defined as the reciprocal of the radius of convergence of the power series z r P r hut 
by a classical result this is precisely p(P) when P is a finite complex matrix, see l fW2| , 



Theorem 5.5], or see the more general Theorem by Hensel |W2, Theorem 5.4]. We 
also call a nonnegative left or right eigenvector corresponding to the Perron root of a 
nonnegative matrix a Perron vector. We usually denote the Perron root of a projection 
matrix P by r and left and right Perron vectors by v 1 and u respectively. 



Since the eigenvalues and eigenvectors (with proper normalization) are continuous in 
the elements of a matrix and every nonnegative matrix is a limit of a sequence of 
irreducible nonnegative matrices, one immediately deduces: 

Corollary 2.2 Let P be a nonnegative matrix. Then 

(a) The spectral radius p(P) of P is an eigenvalue of P with a corresponding nonneg- 
ative Perron vector. 

(b) The spectral radius of the matrix P does not decrease, resp. increase, if any entry 
of it increases, resp. decreases. □ 



An irreducible nonnegative matrix P is primitive if there is only one eigenvalue of 
P that attains the modulus p(P). In this case, the Perron root is properly called 
the dominant eigenvalue of A, which is the usualterm in the demographic literature 
(even under less restrictive conditions), e.g. fCas2, p. 83]. A nonnega tive matrix P 
is primitive if and only if P k is positive for some positive integer k, l |HJ| , Theorem 
8.5.2]. This shows that primitivity (like irreducibility) depends only on the pattern of 
a nonnegative matrix, i.e. if P is primitive (irreducible) then every matrix that has 
positive entries in exactly the same positions is also primitive (irreducible). Applying a 
standard result on themnvergence of powers of matrices l |W2| , Theorem 3.5], one may 
derive from Theore m the following well known result, for closely ^rela ted results see 
for^xample l |Cas2| , p. 86], or in the special case of Leslie matrices, fPal, Lemma 4.5.1] 
or ||, where proofs are given in the same spirit as ours. 

Theorem 2.3 Let P be a primitive nonnegative matrix with spectral radius p(P) = r 
and left and right Perron vectors v t and u such that v t u = 1. Then 



lim (P/rf = uv l 

fc^oo 



(6) | lim 

□ 



The follo wing a pplication of Theorem fc.3| is called the fundamental theorem of demog- 
raphy in IJCush, Theorem 1.1.2], where a different proof is given. 

Theorem 2.4 LetP £je the projection matrix of a standard population model Xk, k = 
0, 1, . . . , given by (Hi). Suppose that P is primitive with spectral radius p{P) = r and 
has left and right Perron vectors v l and u resp. normalized so that v t u = 1. Then 



lim Xj~/r 

k^oo 



[V X U. 



Consequently, if\w\ denotes the sum of entries of the vector w, so that \xk\ will denote 
the total population at time k in the population model, then 



lim \xk\ 

fc— >oo 



'0 if r < 1, 

|(v*a:o)ii| if r = 1, 
oo if r > 1. 



□ 
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Mathematically, a population is a vector that is nonnegative and nonzero. As is usual 
in the demographic literature, we call a population x stable (for a given matrix model), 
though this terminology is not consistent with definitions of stability in other parts of 
mathematics. If, for some positive r, Px = rx, and we call x a stationary population if 
Px = x . We call a population eventually stable if limfc_ +QO Xk/r k exists and is nonzero, 
and we call the population eventually stationary if lim^oo x^ exists and is nonzero. 
The spectral radius p(P) =f p sis called the growth rate of the model, which for primitive 
P is justified by Theorem |2.4| . 

Suppose that P is primitive. It follows immediately from the fundamental theorem 
that, whatever the initial population xq, the number of individuals in x^ grows to 
infinity if r > 1, shrinks to if r < 1, and remains finite if r = 1. Furthermore, in all 
cases, there exists a unique stable population (except for a constant factor), which is 
a stationary population if r = 1. Furthermore, in this case, the fundamental theorem 
shows that every population is eventually stationary. 

The assumption that P is primitive cannot be omitted from this last remark, for the 
conclusion (gj of Theorem |2.4| depends strongly on this assumption. If one merely 
assumes that P is irreducible one may show that Xk/r k , k = 1, 2 . . . is bounded above. 
For reducible nonnegative P the description of the possible limiting behavior of P k XQ is 
quite complicated, see |FS|| for applications of graph theoretic concepts to this problem. 
An example is given near the end of this article where r = 1, but all populations except 
for the stable populations grow to be infinitely large. 

For primitive P, the reciprocal of the growth factor r of P may also be interpreted as 
a factor for stationarity for the model, viz. if both T and F are scaled by the same 
factor 1 jr then the resultant model with matrix P' = (T + F) jr has the property that 
every population is eventually stationary. We are however particularly interested in 
scaling the fertility matrix F without scaling the transition matrix T so that in the 
resultant model every population is eventually stationary, and this leads naturally to 
the considerations in the rest of this paper. 



3 The net reproductive rate 

The main result of this sectionals a somewhat stronger formof a comparison theorem 



due to Cushing and Yicang, IjCY, Theorem3], see also I fCush , r igorem 1.1.3]. Using 

we give a very short 
In the next section 




standard results of Perron- Frobenius theory reviewed in Section, 
and simple proof of this theorem stated below as Theorem 
we show that this theorem can be further refined and generalized using some more 
nonnegative matrix theory. 



Theqrerxi 3.1 Suppose a standard matrix model of population dynamics satisfies (|T 



and 



ul 



and assume 



[hat the projection matrix P 



T + F is irreducible where T is 



nonzero and satisfies (|]j. Denote the growth factor p(P) by r and the net reproductive 
rate p(Q), where Q = F(I — T) _1 , by Rq. Suppose that Rq > 0. Then 



p(T + F/Rq) = 1, 



(7) 



stable 



and one of the following holds: 
r = Rq = 1, or 



1 < r < Rq, or 



< Rq < r < 1. 



(8) compl 



Proof: Since p(T) < 1, by the matrix (I — T) 1 is nonnegative and hence so is 
Q = F(I - T)- 1 . Clearly F^O since p(Q) > 0. 

a rable J EE2 
), note that by Corollary \z.'4 , there exists a nonnegative left 

eigenvector y* of F{I — T) _1 corresponding to the eigenvalue Ro, i.e. y t F(I — T) _1 = 
RoyK Then y l F = Roy\l-T) and hence y t (T+F/R ) =y& Since i? > 0, the matrix 
T + F/Rq is irreducible and hence it follows by Theorem ^.i|(b) that p(T + F/Rq) = 1. 

icpmp 1 

To prove that one of the conditions in (§) holds we consider three cases. 

(i) If R = 1, then 1 = p(T + F) = r. 

(ii) If Ro > 1, then 

T + F/i^o <T + F < R T + F 
with equalities excluded since -F is nonzero. Hence by Theorem 

1 = p(T + F/Rq) <p(T + F)=r< p(R T + F) = R . 



(hi) If < Ro < 1, then again by Theorem fc.l| (c), 

1 = p(T + F/Ro) > p(T + F)=r> p(R T + F) = R . 

□ 

Icga gO 

The special case of Theorem |3.1| when T is strictly lower triangular and F is upper 
iangular is known in numerical linear algebra as the Stein-Rosenberg Theorem, see 



V], pp.68-70 ^for a proof to which our proof o jP, the^ more general result is somewhat 



similar. See 



for a result close to Theorem b.ll. 



rnnripO 



Corollary 3.2 Under the hypotheses of Theorem J7I, consider the modified model 
given by P = T + F/Rq with left and right Perron vectors v t and u such that v t u = 1. 
If P is primitive, then for every initial population xq we have 

lim (T + F/R ) k x = (e*x )«. 

k— >oo 



Proof. This follows immediately by Theorem b.l| (aj and Theorem ]2.d[ , since P is also 
primitive. □ 



There is a corresponding theorem for general nonnegative matrices. 

Theorer^i 3.3 Suppose a standard matrix mo 4^ r ^f population dynamics satisfies 
and (yj , and the transition matrix T satisfies (||). Denote the growth factor p(P) by r 
and the net reproductive rate p(Q), where Q = F(I — T) _1 , by Ro- Then one of the 
following holds: 



ul 



r = R 
If Ro > 0, then 



or 1 < r < Rq, or 



p(T + F/Rq) = 1. 



< Ro < r < 1. 



(9) compgl 

(10) compgO 
□ 



Most of the derivation of this theorem by continuity from Theorem |3.1| is straightfor- 
ward and therefore omitted. However, we shall show why r = 1 implies that Rq = 1. 
For positive e, let F(e) = F + sE, where E is a matrix of the appropriate size all of 
whose entries are 1. Let P (e) = TA-Fj^.). Since P(e) is irreducible, we have p(P(e)) > 
p(P) = 1, and hence by Theorem p7i|, p(F(e)(I — T) _1 ) > 1. If r = 1, letting e tend to 
we obtain R > 1. Now consider P'(e) = P{e) / p{P{2e)). Since p{P{2e)) > p{P{e)) 
we have^^e)) < 1 and hence p{F / p{P{2e))\l - T / p{P{2e)~ 1 ))) < 1, again by The- 
orem 3. ij. Letting e tendto we now obtain Rq < r = 1. It jojlows that Rq = 1. We 
observe that Theorem js^may also be deduced by means of I Sch, Theorem 4.5]. 



4 The refined stability and comparison theorem 

In this section we show that the hypothesis Rq = p(Q) > in Theorem actually 
follows from the remaining assumptions thus allowing us to state a refined version of 
thia^heorem. Our proof depe nds on t^e following proposition, which is a restatement 



officii, Lemma 3.4], see also |Sz] and 



Proposition 4.1 Let T and F be nonnegative matrices with p(T) < 1 and F ^ 0. 
Suppose T+F is irreducible and Q = F(I— T) . Then, after a permutation similarity, 



Q 



Qn Q12 
o o 



(11) splitO 



where Qn is a nontrivial irreducible nonnegative matrix, Qyi is a nonnegative matrix 
every column of which has a positive entry, and the rows of Q correspond to the 
rows of F, if any. □ 



Our stability and comparison theorem may now be stated as: 

Theorern 4.2 Suppose a standard matrix model of population dynamics satisfies (jf^ 



and (gJ7 and^ssume that the projection matrix P = T + F is irreducible with T 
satisfying and F / 0. Denote the growth rate p(P) by^^rjd the net reproductive 
rate p(Q), where Qi^ ~ T)~~^, by Rq. Then equation (/]) holds and so does one of 



the conditions in Furthermore, the matrix Q is irreducible if and only if every row 
of F contains a positive element. 

P-^oofi We have p(Q) = p{Qu) and since Qn is irreducible, p(Qn) > by Theorem 



ill "(which was also observed in fjz|). □ 



T^he first conclusion followsufrom Theorem p.l| and the second from Proposition 



Seneta |Sen| , p. 42] gives a proof of the Stein-Rosenberg theorem which may be adapted 
to show Rq > under the hypotheses of Theorem |4.4 The key ingredient is the 
observation that (/ — T) -1 and (/ — T/r) _1 have the same zero-nonzero pattern. 

In contrast to the situation for irreducible P, one may construct examples of a reducible 
projection matrix P with a nonzero fertility matrix F such that the corresponding net 
reproductive rate Rq = 0. In this connection, we shall prove the following theorem. 

Theorem 4.3 Let P, T and F satisfy the hypotheses of Theorem \j^ V ^Then the net 

reproductive rate Rq > if and only if for some a > 0, p(T + aF) > p(T). 
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Proof: First suppose that Rq > 0. Then, for sufficiently large positive a, we have 
aRi 



o 



> 1. But this is the net reproductive rate of the projection matrix T + aF and 
hence, by Theorem |5.3|, p(T + aF) > 1 > p(T). 



Conversely, suppose p(T + aF) > p(T) for some positive number a. Let det(A7 — 
(T + aF)) = Ylk=i ^ k fk{o). Using the usual determinantal expansions, we see that 
the fk(a), k = 1, . . . , n, are polynomials in a. But they are also signed sums of the 
A;-th elementary symmetric functions of the eigenvalues of the matrix T + aF. Since 
p{T + aF) > p(T), there exists k such that /fc(a) is not a constant polynomial, and 
hence, for this k, \fk(a)\ is unbounded as a goes to infinity. Hence at least one eigen- 
value cannot be bounded in a, and it follows that p(T + aF) is unbounded, and thus 
there exists a such that p(T + aF) > 1. Again applying Theorem T2io T + aF, we 
obtain aRo > 1 and hence Rq > 0. □ 



Our proof of Theorem O~also shows that if Rq is positive, then F may be scaled to 
achieve an arbitrarily large growth rate in the general nonnegative case. Returning to ^ 
an irreducible projection matrix P = T + F with F^Owe observe that Theorem [4.3| 
provides a second proof that Rq > 0, since in this case p(T + aF) > p(T) for all a > 0. 

We now generalize Theorem [op lf F + T is irreducible with Perron root r and positive 

left Perron vector z t then it easily follows that„z*F(J — T/r)^ 1 = rz l and since z* > Q, . 

, *^ HJ— . v ' ' Tcompfin 

wc deduce that p(F(l - T/r)- 1 ) = r, see e.g. If pj| , Cor. 8.1.30]. Thus Theorem |4.2| is 
the special case s = r, q(s) = 1 of the result which now follows. 



Theorem 4.4 Let P, T and F satisfy the conditions of Theorem 
define 



^^or s > p(T) 



Then q(s) > 0. Let P{s) 
net reproductive rate is 



q(s)=p(F(L-T/ S )- 1 )/s. (12) 
T + F/q(s). Then its growth rate, p(P(s)), is s, and its 

R (s) = R /q(s). 
Further, one of the following holds: 

1 = s = Ro(s), or 1 < s < Rq(s), 

Proof. We observe that Q(s^=_F^I 



or < Rq(s) < s < 1. 



T/s)~ l is nonnegative since p(T/s) < 1 and 
hence q(s) > by Theorem Bj^s as F + T/s is irreducible. Let z t be the left Perron 
vector of Q(s). Thus z t Q(s) = z t F(L — T/s)^ 1 = sq(s)z t . An easy computation now 
yields 



z\T + F/q(s)) 



sz 



and hence s is the Perron root of T + F/q(s) as asserted in the theorem. The corre 
sponding net reprodurtjvj^ite is Rq(s) = p((F/q(s))(I — T) _1 ) = Ro/q(s). Then ( |T3| 



assign 



(13) compla 



lpla 



follows by Theorem 



□ 



pla 



We note that the inequalities (|l3|j are strict and are equivalent to p{F(L — T) ) < 
p(F(I - T/s)- 1 ) if p(T) < s < 1 and p(F(I - T)- 1 ) > p{F(L - T/s)" 1 ) if s > 1. 
Further, since p(F(I — T/s) -1 ) is a decreasing function of s, it follows that q(s) is a 
strictly decreasing function of s and that lim^oo q(s) = 0. 



5 Demographic interpretations 



We now discuss the demographic interpretations of our theorems. As we shall impose 
various assumptions on the projection matrix P, we begin each paragraph with an 
assumption which holds throughout the paragraph. 

Assume that P is a primitive nonnegative matrix. As previously observed, the recipro- 
cal of r = p{P) may be viewed as a factor for the scaling of the fertility and transition 



stationary 



matrices in order to obtain a model with every initial population eventual 
population, viz. p(P/r) = 1. Similarly, since Rq = p(Q) > 0, (Theorem Of/" we may 
leave the transition matrix fixed and scale the fertility matrix by the reciprocal of Rq 
and in the resultant model every initial population will tend to a stationary population 
which is unique except for a multiplicative constant (Corollary \i.'4 )- Intuitively, scal- 
ing the fertility matrix only should be more radical than scaling both the fertility and 
the transit jion mat rixt o achieve the same objective, and that is exactly the content of 
Theorems b.l and 



Suppose now that P is irreducible. It is possible to adapt the previous paragraph to 
this case by considering scalings to achieve a growth rate of 1, but we shall immediately 
turn to interpeting the more general Theorem pf.4| . Given any s > p(T), it is possible 
to obtain a growth rate of s without changing the transition matrix T by scaling 
the fertility matrix F to F/q(s), where q(s) is given by (\l'4). since q(s) is a strictly 
decreasing function of s for s > p(T), to achieve a higher growth rate we require greater 
fertility, which again is intuitively clear. We shall discuss these points further when we 
turn to the Leslie model. 

Suppose that P is irreducible. By Proposition [Sj every column of the submatrix 
(Qiij Q12) has a positive entry, and this implies that every class in the injtia| population 
has descendants during its lifetime. By the last assertion of Theorem ^.4 the matrix Q 
is irreducible if and only if there are newborns in every population class in the case of 
a population that has members in each class. 

Now assume again that P is irreducible. We turn to characterizations of Rq which may 
be obtained from Theorem |4.2[ and from the following characterization of the Perron 
root of an irreducible nonnegativernatrix A which is a reformulation of a well-known 
characterization due to Wielandt IflWlfl, see also 



p.65]: 

p(A) = max{s : Ax > sx, for some x > 0, i / 0} 
= min{s : Ax < sx, for some x > 0, 1^ 0}, 

and the equality is attained in either inequality if and only if x is the Perron vector 
of A. For this formulation and references to it, see the commentary following |W1 | 
in Wielandt's Mathematical Works. Let z denote a nonnegative vector with the same 
number of entries as there are columns in Qn, and let Z{ be its ith entry. Then there 
exist indices i and j such that (Qnz)i > R§Zi and (Qnz)j < R^Zj. We may put the 
case r = Rq = 1 into words thus: For every initial population, there must be one 
class of newborns that over its lifetime produces at least as many descendants in the 
same class and there must be one class of newborns that over its lifetime produces at 
most as many descendants in the same class. Furthermore, there exists a population 
of newborns (i.e. corresponding to the nonzero rows of F) which over the course of its 
lifetime reproduces itself exactly if and only if there exists a stationary population (but 
note that this population of newborns is in general not part of a stationary population). 



Q 



cible, it is easily pr. 
H Theorem 2.20], 



fed using an additional 
Corollary 8.4.6 and 



Still under the assumption that P is irr^ 
part of the Perrpn-Frobenius theorem, 
Remark 8.4.9], |V|, Theorem 2.3], l|G], Theorem 2, p. 53], that for some positive integer 
d, P d is the direct sum of primitive matrices, and hence that there is a population that 
is stable in a periodic sense for this model, that is for some positive integer d and every 
initial population xo there exist populations wq, . . . , Wd-\ (depending on xo) such that 
(Hindoo Xkd+i)/ rk = w ii i = 0, ■ ■ ■ ,d — 1. To obtain the previous primitive case, we 
put d = 1. Another equivalent condition for p{P) = 1 in the irreducible case is that 
there should exist an initial population xq ijy^£h that lim^oo xt exists and is nonzero. 
This may be proved by means of Theorem 

We now turn to the case of the Leslie model, where T is a matrix with nonzero elements 
on the first subdiagonal and O' s^ els ewhere, and F is., a matrix all of whose nonzero 
elements are in its first row, see 




s^lsewrgre, and t i s^ a m atrix 
fCesl| ], |Pa| p.38] or | |Cush| , p.4]: 



h 












and 



F 



tn-l 0/ 



(h 




fn-l 




fn\ 




0/ 



Whether P is 
newborns 



sec 



^re ducible or not, the situation is simpler as there is only one class of 
Par] for a discussion of reducible Leslie models. In this case, the net 
reproductive rate Ro is the (1,1) e ntry o f Q and equals the sum of the elements of 
the first row of F, as observed in I flCush , p. 9, 1.7] where it is remarked that the net 
reproductive rate is "the expected number of offspring per newborn over the course of 
its lifetime". Thus in this case the no np ro ducti ve ratagquals the net reproductive 



rate as defined in the early papers l |Bern |, f Cew ] and f Ces2j , p. 234]. In turn, this 



an adaptation of a conc ept pr eviously used in continuous population models, see 
jou, p. 11 5]. By Theorem 373 we have Rq = 1 if and only if r = 1, and if P is 
a primitive projection matrix for a Leslie model then this has the interpretation that 
every population is eventually stationary if and only if the expected number of offspring 
over the course of a newborn's lifetime is 1, viz. the (1, 1) entry of Q is 1. 

Let P = T+F again be a Leslie matrix. Then q(s) is the leading entry of F{I—T/s)~ l /s 
which is 



q{s) = /is" 1 + f 2 tis- 2 + ■■■ + /„(tn_i • • • ti)s" n , 

a polynomial in s -1 with coefficients involving all nonzero entries in F and T. 
classical that the nonzero eigenvalues A of P satisfy the equation q(X) = 1, e. 



(14) I poly 




p. 42], and this is consistent with our more general results since by Theorem O the 
growth rate r of P satisfies q{r) = 1. We note that the operation of s caling the fertility 
matrix by a constant factor was previously considered by Leslie in f L~es2| , Sec. 5(b)]. 
His results in this area have a somewhat complicated appearance since he wished to 
determine s as a function of q(s), which requires the solution of a polynomial equation. 
This gives rise to the following observation: To determine the growth rate for a Leslie 
matrix one needs to find a positive root of a polynomial equation, but to scale the 
fertility matrix to achieve an assigned growth rate s it suffices to divide the fertility 
matrix by q(s) in (|i4|j. 

Care needs to be taken to interpret the various terms used in this article when the 



in 



matrix P is reducible. It is instructive to consider a simple example. Let 



P 



1 
1 1 



For the projection matrix P, an initial population is stable if and only if its first 
element is 0, otherwise it will tend to infinity. The growth rate r = p(P) = 1, and the 
net reproductive rate Rq = 1 whatever ^avb^ the Jra^isition matrix T and fertility 



matrix i 7 chosen subject to conditions (|1|), (^[) and 



We conclude by giving a numerical exampJewhich is based on a plant lifecycle involving 
vegetative as well as seed reproduction, IjCasl, Example l.c]. Let 



;i/2) 



/0 0\ 

1 

1 

110 

\0 1 0/ 



and 



F = (1/2) 



/0 1\ 



1 



\o 0/ 



Let P = F + T. This matrix is irreducible, but imprimitive. Then the growth rate 
r = V2/2 and the stable populations for P are u = (a/2, 1, 3, 2\/2, 2)* and its positive 



multiples. Further the next generation matrix is 



Q = (1/8) 



/l 1 1 2 4\ 



2 2 2 4 



\0 0/ 



indicating that the population of newborns that reproduces itself with the same distri- 
bution of newborns in the next generation is w = (1,0,2,0,0)*, the Perron vector of the 
next generation matrix Q. Note that the submatrix in rows and columns 1 and 3 of Q is 
irreducible and that R = p(Q) = 3/8. We also note that (P/p(P)) k w, k = 0, 1,2, . . . 
does not tend to a limit, which indicates that generational stability is compatible with 
permanent oscillation of the normalized population distribution over time. 

If P is replaced by Pi = T + 8F/3, then correspondingly Q x = (8F/3)(I -T)' 1 = 
and we have both p{P\) = p(Qi) = 1- Since p(T) = 0, given any s > 0, we can use (|T2| 
to compute q(s) such that P(s) = T + F/q(s) has growth rate s. We obtain q(s) = (1 + 
2s 2 )/8s 4 , and hence the corresponding net reproductive rate is Ro(s) = 3s 4 /(l + 2s 2 ). 
The corresponding stable population is (4s 3 , 2s 2 , 2s 2 + 8s 4 , 2s + 4s 3 , 1 + 2s 2 )*. 
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